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Abstract 

We study the creation of disoriented chiral condensates with some initial 
boundary conditions that may be expected in the relativistic heavy ion col- 
lisions. The equations of motion in the linear a- model are solved numerically 
with and without a Lorentz-boost invariance. We suggest that a distinct cluster 
structure of coherent pion production in the rapidity distribution may emerge 
due to a quench and may be observed in experiments. 
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Recently much attention has been paid to an interesting proposal to observe a 
disoriented chiral condensate (DCC) in high energy collisions fj], @, ||, ||. Such events 
can be signaled by a coherent pion emission along some particular isospin direction 
in the collision domain. A probability distribution of neutral pions is characterized 
by a non-binomial function if all isospin directions are equally probable, 

pm = (i) 

where r = n n o / (n^o + n n ± ) . This behavior may have already been observed in the 
so-called Centauro events in cosmic ray collisions ||. The central theoretical work is 
to understand the mechanism for creating such a DCC domain in relativistic heavy 
ion collisions. Rajagopal and Wilczek || have suggested that the expansion of highly 
relativistic debris from collisions "quenches" the high temperature field configurations 
such that the long wavelength modes will grow with time, leading to a correlated 
domain. It is recently reported, however, that the correlation size may be very small 
0- 

In this paper, we shall study the creation of the disoriented chiral condensate with 
a more realistic initial condition that may be expected in the high energy collisions. 
We shall re-examine the idea of a quench based on a linear a- mo del and suggest that 
indeed, following a quench, correlated pion fields will populate near the light cone 
where the leading collision particles are expanding at speed of light. The production 
of the unusually rich (or poor) 7r°'s will cluster into groups in the rapidity distribution. 
We are inspired by the consistency of our results with some recent events of JACEE 
experiments || and confirm the general picture of a "Baked Alaska" sketched by 
Bjorken, Kowalski and Taylor ||. 

We assume that following a quench the possible DCC is to be described by classical 
low energy effective interactions of pions at zero temperature 

£ = J d*x S^d^d^i - j (0,0; - v 2 ) 2 + Fa | (2) 

where 0j = (er, 7r) stands for a vector in internal space. Ha is an explicit chiral 
symmetry-breaking term which is responsible for the mass of the pseudo-goldstone 
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bosons, the pions. We look for a solution to the classical equations of motion of 
with a given initial condition. Clearly the initial condition is essential to determine 
whether or not there will be DCC in the system. For example, if the and initially 
align in the a direction throughout the space, then the system does not have any 
pion field at any time. What might be the typical initial configurations? In a quench, 
the system begins at a temperature well above the chiral phase transition point T c . 
There are thermal fluctuations in all internal directions. Therefore, it is appropriate 
to assume that initially the space averages ((ft) ~ and ((f)) ~ but (0 2 ) ^ and 
((f) 2 ) both for a and tc fields. In addition, we are modelling a situation applied 
to highly relativistic collisions. The system should satisfy the space-time geometry 
of such a collision where the incident nuclei are, shortly after the collision, highly 
Lorentz-contracted "pancakes" receding in opposite longitudinal direction from the 



collision point [10]. An approximate 1 + 1 Lorentz invariance of the system is indeed 
inspired by the existence of a cental-plateau structure in the rapidity distribution of 



produced particles in high energy cosmic-ray events and pp or pp collisions |T2|]. Of 
course, the above picture should be modified at large transverse distance, comparable 
to the nuclei radii. However, in most of the transverse area, the field changes rather 
slowly in the transverse directions compared with the rapid longitudinal expansion. 
Bearing this approximation in mind, we shall assume a uniform field in the transverse 
direction, i.e. is only a function of t and z. In terms of more convenient coordinates, 
the proper time r = y/t 2 — z 2 and the (spatial)rapidity variable 77 
equations of motion read, 
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Given an initial condition, Eqs. (^) and @ 
evolution of the system. 

The simplest case to solve Eqs. (H) and 
Lorentz-boost invariance, i.e., 



-Act (a 2 + tt 2 - v 2 ) + H ; (3) 

_ A7r (a 2 + tt 2 - v 2 ) . (4) 
completely determine the space-time 

@) is when the system has an initial 



(5) 
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where 0o is independent of rj. The lack of in the Klein-Gordon operator guar- 
antees that Eq. (|5|) is sufficient to maintain the boost invariance at any time (there 
is no need to impose ^r0(r o ,^) = constant). To model a quench, we assume that the 
system is initially in a symmetric phase and lies on the top of the "Mexican Hat" of 
the potential, 

7t(t ) = , <7(t )~0. (6) 

To obtain a non-trivial solution for ir field, we allow a small initial "kick" of the 
system, i.e. a non-zero velocity §^(to) 7^ 0. In this case, since the system starts from 
one point z ~ at t ~ (the rapidity dependence of decouples and the motion is 
zero-dimensional), the correlation in isospin direction is automatically 100%. That is, 
if the initial kick is in 7r°-direction, then = at all time such that there will be 
only 7T° production. A numerical solution is plotted in Fig. 1 where we have used the 
following standard parameter input: v = 87.4 MeV, H = (119 MeV) 3 and A = 19.97 
so that f n = 92.5 MeV, m n = 135 MeV and m a = 600 MeV. The initial condition 
we used in Fig. 1 is 0(r o ) = 0, <90(r o )/<9r = (1, 5, 0, 0) MeV/fm, at r = 1 fm/c. The 
general feature of the solution is following. The a field grows from zero and takes 
about 1 fm/c proper time to reach the true vacuum expectation value (a) ~ f n , while 
the 7r field oscillates around zero rather slowly and eventually tends to zero when the 
proper time gets large r = r max ~ 30 fm/c. In terms of the real space-time variables t 
and z, the disoriented chiral condensate develops in the interior region of the receding 
pancakes, surrounded by the normal vacuum in the exterior; as time (t) evolves, the 
DCC will be squeezed toward the light cone, leaving behind a normal vacuum where 
(tt) = in the central region inside the light cone. For t ^> r max , the it field is found 

2 

only in a small region Az ~ near the light front while the a field occupies most 
of the space inside the light cone, surrounded by a thin shell of DCC. This picture 
supports a space-time geometry of DCC sketched in |J. 

While the essential feature of the above solution may be generic, the real situation 
is more complicated. One complication is that the system may not start from one 
point shortly after the collision and there may not be a Lorentz-boost invariance in the 
system. In this case, the initial condition like Eq. (||) is not appropriate. If the system 
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initially has a distribution in rapidities and is not correlated in isospin directions, the 
the crucial question would be whether a quench can yield a coherent pion field at a 
later time. 

We assume that shortly after the collision, the typical configuration is that of a 
thermal random fluctuation at high temperature. The fluctuation is governed by the 
temperature. Similarly to Ref. H, we choose 4>(to) and <90(to)/cV randomly according 
to gaussian distributions so that 

(^o,*7)>,= <^(t 0) »7)>, = (7) 

but (<f) 2 (T ,i])) v = v 2 /A and ((^0(r O) r])) 2 ) v = v 2 . By solving equations (|) and (|), 
we find that as the proper time evolves, the average value (a) v grows from zero quickly 
and takes its final value (unlike the case with boost-invariant initial condition, the 
oscillation of (a) v around /„- is quickly damped at large r); (ir) v grows from zero and 
oscillates around zero slowly in the same way as in Fig. 1. Again, (71")^ decreases to 
zero at about r ~ 20 ~ 30 fm/c. This clearly indicates a strong correlation among 
the isospin directions in rapidities. Indeed, we can define a correlation function at a 
given r, 

cm-,o= ; (T -?r ( ;;"l ■ «» 

TT Z (T, 1]) + 7T Z {T, 1]') 

We find that after about 1 fm/c evolution in proper time (which is the time scale 
that it takes for the average value of a field to reach its final value {a) v ~ /„-), 
C(t,i] — rf) changes from a zero value corresponding to a typical random distribution 
to an exponential distribution with non- vanishing width as shown in Fig. 2. For small 
At] — Tj — 7/, C(t, Arj) can be fit by an exponential function 

C(r,Ar]) oc exp(— c m n T\ Ar]\) (r — t > 1 fm/c) (9) 

where Cq ~ 0.4. At r > 20 ~ 30 fm/c, the correlation disappears as indicated by 
the average value (71%. A similar result is derived in based on a non-linear o- 
model in 1 + 1 dimensions where the pion fields are the phases of the order parameter. 
However, in our model, we find that if initially (<f)) 2 ~ f 2 , a correlated distribution 
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does not occur. This clearly shows that the emergence of a correlation in isospin 
orientation is mainly due to the unstable long wavelength modes of the pion fields when 
(4>) 2 < f%. Our result unambiguously points to the success of the "quench" mechanism 
for creating DCC (i.e. a high temperature initial condition plus zero temperature 
equations of motion). 

Eq. (H) indicates that a rapidity interval in which the pion field is correlated in 
isospin direction could be as large as 2 ~ 3. Our conclusion coincides with the recent 
cosmic ray experiment || where the dominant neutral pions are found in an interval 
A 77 ~ 3. Eq. (||) also suggests that the correlation occurs mainly in the region near 
the light cone where r is small. At a given large t, a region of ordinary vacuum where 
(77)^ vanishes effectively is found inside a shell of coherent pion field. However, the 
cluster structure of pions radiated from the coherent pion field may occur anywhere 
in the whole rapidity region. We should also emphasize that the width of the cluster 
size in our study depends on the initial time Tq when the quenching happens. For 
sufficiently large initial time, the width can become very small. 

Although the size of correlated DCC domains may be small in space, especially at a 
large t because of the contracting effect due to the relativistic expansion, we conclude 
that the DCC production can be distinctly observed as a cluster structure in rapidities. 
What is less clear is whether there is indeed a quench after the collision. Strictly 
speaking, the quench condition is an idealization of the more complex situation where 
the temperature relaxes to zero only gradually. For example, the hydrodynamics 
suggests a temperature drop according to T = To(^) 1 / 3 flQfl . In this case, equations 
(|3|) and may not be appropriate. Clearly, more theoretical work is needed to 
examine if DCC can be created in the process of cooling. 
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Figure Captions 



FIG. 1. Proper time evolution of a and tti fields following a Lorentz-boost invariant 
initial condition at r = 1 fm/c. 

FIG. 2. The correlation function C(r, A77) of DCC. The initial configuration is 
that of a thermal fluctuation at r = 1 fm/c. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312005vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312005vl 




Fig. 2 



